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ABSTRACT 

Aims. We study the variability of the nova-like cataclysmic variable TT Ari, on time-scales of between minutes and 
months. 

Methods. The observations in the filter R were obtained at the 40-cm telescope of the Chungbuk National University 
(Korea), 51 observational runs cover 226 hours. The table of individual observations is available electronically. For the 
analysis, we have used several methods: periodogram, wavelet and scalegram analysis. 

Results. TT Ari was in the "negative superhump" state after its return from the "positive superhump" state, which lasted 
8 years. The ephemeris for 12 best pronounced minima is r„i„ = B JD2453747.0700(47) -|-0.132322(53)£. The phases of 
minima may reach 0.2, indicating non-eclipse nature of these minima. The quasi-periodic oscillations (QPO) are present 
with a mean "period" of 21.6 min and mean semi-amplitude of 36 mmag. This value is consistent with the range 15-25 
minutes reported for previous "negative superhump" states and does not support the hypothesis of secular decrease 
of the QPO period. Either the period, or the semi- amplitude show significant night-to-night variations. According to 
the position at the two-parameter diagrams (i.e. diagrams of pairs of parameters: time, mean brightness of the system, 
brightness of the source of QPO, amplitude and time scale of the QPOs), the interval of observations was splitted into 5 
parts, showing different characteristics: 1) the "pre-outburst" stage, 2) the "rise to outburst", 3) "top of the outbursts", 
4) "post-outburst QPO" state, 5) "slow brightening" . The brightness of the QPO source was significantly larger during 
the 10-day outburst, then for the preceding interval. However, after the outburst, the large brightness of the QPO source 
still existed for around 30 days, producing the stage "4". The corresponding diagram mqpoijn) has a splitting into 
two groups at the brightness range 10^6-10'"8, which correspond to larger and smaller amplitudes of the QPO. For the 
group "5" only, statistically significant correlations were found, for which, with increasing mean brightness, also increase 
the period, amplitude and brightness of the source of QPOs. The mean brightness at the "negative superhump state" 
varies within 10?^3-11™2, so the system is brighter than at the "positive superhump" (11™3), therefore the "negative 
superhump" phenomenon may be interpreted by a larger accretion rate. The system is an excellent laboratory to study 
processes resulting in variations at time-scales from seconds to decades and needs further monitoring at various states 
of activity. 

Key words, novae: cataclysmic variables - stars: variables: general - binaries: general 



1. Introduction 

TT Ari is one of the brightest and thus extensively stud- 
ied nova-hke cataclysmic variable (see e.g. monographs: 
Warner (1995), Hellier (2001)). A wide variety of processes 
takes place in this object, making it variable at time scales 
from 9.6 seconds (Kozhevnikov, 1986) to years (Hudec et 
al. 1984, Bianchini 1990, Kraicheva et al. 1999). Due to a 
low inclination, there are no prominent eclipses at the light 
curve, but there was a "negative superhump" -type mod- 
ulation and ~ 20— min quasi-periodic oscillations (QPO). 
Mardirossian et al. (1980) reported on ~ 40** "oscillatory 
events" , > 6-min QPOs, flickering and a nearly sinusoidal 
shape of 3— hour modulations. 

Krautter et al. (1981) classified TT Ari as "a dwarf nova 
at a permanent outburst". Hoard (2007) included TT Ari 

* Tables 1,2,4 and Figures 7,10 are only available in electronic 
form at ,http://www.edpsciences.org| 



in his "Big hst" of the SW Sex stars (see e.g. Hellier 2000 
for a review). Semeniuk et al. (1987) suggested 4 periods 
in TT Ari, namely, 3.8 day beat period of the 3.2 - hour 
photometric and 3.3 hour spectroscopic periods, and quasi- 
periodic oscillations (QPO) with an apparent decrease of 
the QPO period from 27 (in 1961) to 17 minutes (in 1985). 
Kozhevnikov (1986) found variability at the short timescale 
9.6 sec. 

To study these types of variability, several international 
observational campaigns have been organized, the results 
of which have been summarized by Wenzel et al. (1986), 
Tremko et al. (1996), Skillman et al. (1998) and Andronov 
et al. (1999). Hollander and van Paradijs (1992) suggested 
a dependence of the characteristics of QPOs on accretion 
rate. Kraicheva et al. (1999) suggested variations of the 
accretion rate as a cause of switches between the positive 
and negative superhumps. 
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Fig. 1. The complete light curve of TT Ari for 51 nights of observations (bottom) and of the comparison star "d". The 
V magnitudes in one night (53683) of observations were shifted using the mean value of the instrumental color index 
V-R=a°223. 



Table 2. Table of R observations (HJD-2400000, magni- 
tude) obtained at the 40cm telescope of the CBNU obser- 
vatory, Korea. The table (14199 lines) is published electron- 
ically only. 



One of the most exciting events in TT Ari was a switch- 
ing in 1997 from the state of "negative" superhumps, which 
lasted nearly 3 decades of optical studies, to a state of "pos- 
itive" superhumps (Skillman et al. 1998), in an excellent 
agreement with the "Por5 — -Ps/i" ("orbital period - super- 
hump period") statistical relation (Andronov et al. 1999). 
The system showed highly asymmetric (M-m=0. 320(5)) 
permanent superhumps with a same amplitude of 0™19 
both in V and R and a period of 0^14841(1) in 2004 
(Andronov et al. 2005). However, in 2005, the system has 
changed its state back to negative superhumps with possi- 
bly some additional evidence for instability (Andronov et 
al. 2005). It caused to start a regular photometric mon- 
itoring at the Chungbuk National University. Preliminary 
results were briefly published by Yoon et al. (2006). Results 
of this observational campaign are presented in this paper. 



2. Observations and Comparison Stars 

2.1. Observations 

The time-series observations in R band were made using a 
512x512 FLI CM-1 CCD attached to the 40 cm telescope 
at Chungbuk National University (MAEDE LX200). The 
field of view of a CCD image is about 12. 97x 12.97 arcmin^, 
given a CCD plate scale of l."52 per pixel at the f/6.3 
Schmidt-Cassegrain focus of the telescope. 

Instrumental signatures of each CCD frame were re- 
moved and calibrated using the bias, dark, and flat field 
frames, with the aid of the IRAF package CCDRED. We 
obtained instrumental magnitudes of stars from the empiri- 
cal point-spread function (PSF) fitting method in the IRAF 
package DAOPHOT (Stetson 1987; Massey & Davis 1992). 

The journal of observations is presented in Table 1. All 
Julian dates are expressed omitting " 24" at the beginning, 
so the integer part contains 5 digits instead of 7. For the 
legend of a separate night of observations, we have used the 
integer "legend Julian Date" LJD=int(BJDi — 0.5), where 
BJD is bary-centric Julian date, BJDi is BJD for the first 
observation, and int(a;) is an integer part of x. An additional 
substraction of 0.5 was needed because the integer part of 
BJD changed during the night many times. This method 
prevents duplication of LJD for different runs. E.g. the run 
starting on 53684.04 will have a designation LJD=53683, 
whereas the next run starting at 53684.97 (LJD=53684). So 
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Table 1. (Available electronically only). Journal of observations of TT Ari: The "legend Julian date LJD" (integer part 
of the Julian date for the beginning of the night); number of observations n; time of the begin tb and end tg of observations 
(for the many nights, the observations ended on the next integer JD); duration of the run te — h; magnitude range for 

individual data points m„iax, "nimin] nightly mean m and it's accuracy estimate a{fri): r.m.s. deviation of the single 
observation from the mean a{m)] cxp - time resolution in seconds. All moments of time are expressed in (BJD-2400000). 



cr(m) a(rn) exp 



LJD 



n 



tb- 



duration 



range 



m 



53668 
53682 
53683 
53684 
53685 
53686 
53691 
53692 
53693 
53694 
53696 
53698 
53700 
53703 
53704 
53706 
53707 
53712 
53714 
53718 
53719 
53721 
53723 
53727 
53728 
53729 
53730 
53731 
53733 
53741 
53742 
53743 
53753 
53756 
53757 
53758 
53759 
53760 
53761 
53762 
53767 
53769 
53774 
53776 
53777 
53778 
53783 
53790 
53792 
53797 
53800 



89 
127 
243 

213 
326 
445 
545 
371 
302 
331 
326 
376 
38 
159 



391 
31 
301 
408 
53 
60 
463 
419 
325 
138 
418 
467 



426 

490 
302 

355 
201 
428 
403 
403 



228 
177 
144 
289 



53669. 
53683. 
53684. 

53684, 
53685, 
53686, 
53691, 
53692, 
53693. 
53694. 
53696. 
53698, 
53701, 
53704, 



216 53705 



53707 
53707 
53713 
53714 
53718 
53720 
53721 
53723 
53728 
53728 
53729 
53730 



457 53731 
199 53734 



53741 
53742 
53743 
53753 
53756 
53757 
53758 
53759 



87 53760 

199 53761 

305 53762 

217 53767 



53769 
53774 
53776 
53777 



357 53778 
298 53783 



218 
173 
91 
171 



53790 
53792 
53797 
53800 



12931-69. 
16957-83. 
03965-84. 

97582-85. 
98716-86. 
91790-87. 
89585-92. 
94350-93. 
97067-94. 
98199-95. 
92302-97. 
88323-99. 
05339-01. 
14776-04. 
07537-05. 
03617-07. 
89166-07. 
00121-13. 
96644-15. 
93715-18. 
16045-20. 
90647-22. 
89227-24. 
01029-28. 
90736-29. 
88847-30. 
88198-31. 
87834-32. 
00208-34. 
93878-42. 
88202-43. 
88083-44. 
88591-54. 
97976-57. 
89002-58. 
89195-59. 
90846-60. 
89155-60. 
89168-61. 
92413-63. 
94171-68. 
90351-70. 
92833-75. 
90179-76. 
90157-78. 
91810-79. 
92239-84. 
91312-91. 
91504-93. 
95336-98. 
91434-00. 



20075 
29230 
24841 

14497 
29891 
27436 
26499 
25535 
27612 
26401 
16962 
29020 
09047 
27062 
24029 
24017 
92107 
16906 
24025 
96768 
20221 
18562 
20482 
20304 
15139 
15052 
16156 
14991 
14173 
34369 
12333 
04232 
07058 
10862 
07701 
07778 
08968 
92635 
98891 
08970 
05934 
05308 
04274 
97372 
05429 
05101 
03073 
02148 
02008 
00170 
98895 



0.07144 
0.12273 
0.20876 

0.16915 
0.31175 
0.35646 
0.36914 
0.31185 
0.30545 
0.28202 
0.24660 
0.40697 
0.03708 
0.12286 
0.16492 
0.20400 
0.02941 
0.16785 
0.27381 
0.03053 
0.04176 
0.27915 
0.31255 
0.19275 
0.24403 
0.26205 
0.27958 
0.27157 
0.13965 
0.40491 
0.24131 
0.16149 
0.18467 
0.12886 
0.18699 
0.18583 
0.18122 
0.03480 
0.09723 
0.16557 
0.11763 
0.14957 
0.11441 
0.07193 
0.15272 
0.13291 
0.10834 
0.10836 
0.10504 
0.04834 
0.07461 



10.847-11.167 
10.995-11.338 
11.045-11.289 

11.000-11.282 
10.997-11.288 
10.929-11.303 
10.986-11.256 
10.975-11.239 
10.982-11.268 
10.972-11.224 
10.962-11.252 
10.773-11.228 
10.956-11.180 
10.804-11.052 
10.712-10.959 
10.582-10.905 
10.660-10.885 
10.265-10.615 
10.071-10.478 
10.233-10.569 
10.251-10.509 
10.150-10.514 
10.510-10.854 
10.407-10.959 
10.597-11.003 
10.534-10.939 

10.596- 11.038 
10.422-10.938 
10.615-11.013 
10.400-10.887 
10.378-10.933 

10.525- 11.000 
10.482-10.827 

10.597- 10.861 
10.582-10.836 
10.557-10.792 
10.450-10.748 
10.507-10.646 
10.540-10.846 
10.618-10.914 
10.552-10.744 
10.606-10.778 
10.548-10.748 
10.529-10.783 
10.523-10.770 

10.526- 10.748 
10.289-10.590 
10.484-10.709 
10.471-10.673 
10.381-10.563 
10.513-10.709 



10.9975 
11.1820 
11.1667 
11.1341 
11.1767 
11.1435 
11.1194 
11.1228 
11.1418 
11.1039 
11.1090 
11.0404 
11.0932 
10.9180 
10.8413 
10.7678 
10.7859 
10.4350 
10.2661 
10.4005 
10.3689 
10.3352 
10.6845 
10.7631 
10.8262 
10.7779 
10.8290 
10.6882 
10.7913 
10.6814 
10.6474 
10.7455 
10.6851 
10.7513 
10.7254 
10.6759 
10.5982 
10.5905 
10.6917 
10.7826 
10.6707 
10.7000 
10.6563 
10.6487 
10.6635 
10.6468 
10.4454 
10.5921 
10.6038 
10.4842 
10.6096 



0.0073 
0.0057 
0.0033 

0.0042 
0.0028 
0.0031 
0.0021 
0.0026 
0.0029 
0.0028 
0.0028 
0.0028 
0.0075 
0.0044 
0.0035 
0.0031 
0.0106 
0.0040 
0.0036 
0.0113 
0.0082 
0.0031 
0.0031 
0.0059 
0.0065 
0.0038 
0.0039 
0.0045 
0.0053 
0.0043 
0.0042 
0.0055 
0.0035 
0.0042 
0.0022 
0.0023 
0.0033 
0.0029 
0.0041 
0.0033 
0.0027 
0.0025 
0.0030 
0.0049 
0.0029 
0.0024 
0.0036 
0.0041 
0.0025 
0.0045 
0.0032 
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0.0619 
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20 



BJD-LJD is elapsed time in days since the legend Julian 
date, and may exceed unity. For all observational runs, 
BJD LJD ranged from 0.88 to 1.34. Totally 14199 obser- 
vations were obtained during 226 hours in 51 nights (JD 
2453669-2453800). 

The original observations (BJD, R magnitude) are pre- 
sented in Table 2 (electronically only). For two nights, there 
was CCD V photometry (one night a series of V, another 



night sequences VRVR... with alternatively changing fil- 
ters). 

2.2. Comparison stars 

The photometric UBV standards in the field have been 

published by Goetz (1985) by linking to the star "c" (= 
GSC 1207 1535) of Shafter et al. (1985). The chart based 
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Fig. 2. The individual light curves of TT Ari. Abscissa is expressed as BJD-LJD, virhcre the "legend JuUan date" LJD 
(-2400000) is shown near each curve. As the integer part of the Julian date JD often changed during our observations, 
we preferred to use the value corresponding to the beginning of the night. The length of horizontal lines between the first 
and second light curves is equal to the orbital period. 
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Table 3. Brightness of the comparison and chock stars ("c" 
and "d", Goetz 1985) compiled from different publications. 



u 


B 


V 


R 


I 




Ref. 




comparison star "c' 


GSC 1207 1535 




12.04 


11.72 


11.01 


10.41 


9.97 


Efimov ct al. 


(1998) 


11.91 


11.68 


10.99 






Shafter et al. 


(1985) 












Goetz 


(1985) 




11.74 


11.05 


10.66 


10.3 


Henden 


(2007) 




11.82 


11.17 






VSNET 


(2005) 






11.1 






AAVSO 


(2007) 






check star 


"d" = 


GSC 1207 1562 




13.38 


12.17 


11.02 






Goetz 


(1985) 




12.20 


11.05 


10.51 


9.97 


Henden 


(2007) 




11.96 


10.86 






VSNET 


(2005) 






11.1 






AAVSO 


(2007) 



Table 4. (Available electronically only). Characteristics of 
the highest peaks at the "wavelet periodogram" S{P) for 
the individual nights: i - mean time of the run, m - sam- 
ple mean magnitude (as in the Table 1), rnqpo - stellar 
magnitude of QPOs, IgP, Igr, S{P). 
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.94 


-1 


.3872 


0, 


,2839 


53707, 


,90587 


10, 


,7859 


12, 


,9399 


-1 


.98 


-1 


.1273 


0, 


,7236 


53713, 


,08494 


10, 


,4350 


12, 


,8533 


-1 


.82 


-1 


.2328 


0, 


,3574 


53715, 


,10349 


10, 


,2661 


13, 


,0893 


-1 


.92 


-1 


.3947 


0, 


,1543 


53718, 


,95241 


10, 


,4005 


12, 


,4503 


-2 


.00 


-1 


.0857 


0, 


,5146 


53720, 


,18131 


10, 


,3689 


13, 


,3453 


-1 


.98 


-1 


.4559 


0, 


,1677 


53722, 


,05972 


10, 


,3352 


13, 


,4575 


-2 


.13 


-1 


.5143 


0, 


,1299 


53724, 


,04834 


10, 


,6845 


13, 


,7616 


-1 


.75 


-1 


.4962 


0, 


,1488 


53728, 


,10643 


10, 


,7631 


13, 


,1026 


-1 


.81 


-1 


.2013 


0, 


,2479 


53729, 


,01869 


10, 


,8262 


13, 


,1744 


-1 


.80 


-1 


.2048 


0, 


,3683 


53730, 


,01934 


10, 


,7779 


13, 


,2913 


-1 


.87 


-1 


.2708 


0, 


,2663 


53731, 


,02173 


10, 


,8290 


13, 


,1531 


-1 


.87 


-1 


.1952 


0, 


,3237 


53732, 


,01320 


10, 


,6882 


12, 


,8525 


-1 


.83 


-1 


,1314 


0, 


,3356 


53734, 


,07165 


10, 


,7913 


13, 


,5340 


-1 


.90 


-1 


.3625 


0, 


,1931 


53742, 


,05811 


10, 


,6814 


13, 


,2361 


-1 


.81 


-1 


.2874 


0, 


,2354 


53743, 


,00344 


10, 


,6474 


12, 


,9647 


-1 


.81 


-1 


.1925 


0, 


,3037 


53743, 


,96000 


10, 


,7455 


13, 


,2835 


-1 


.90 


-1 


.2807 


0, 


,2038 


53753, 


,97373 


10, 


,6851 


13, 


,7188 


-2 


.02 


-1 


.4789 


0, 


,1812 


53757, 


,04199 


10, 


,7513 


14, 


,3629 


-2 


.06 


-1 


.7100 


0, 


,1480 


53757, 


,98316 


10, 


,7254 


14, 


,3314 


-2 


.08 


-1 


.7077 


0, 


,1520 


53758, 


,96880 


10, 


,6759 


14, 


,3272 


-2 


.07 


-1 


.7258 


0, 


,1330 


53759, 


,98530 


10, 


,5982 


13, 


,9531 


-1 


.96 


-1 


.6073 


0, 


,1150 


53760, 


,90784 


10, 


,5905 


14, 


,5697 


-2 


.36 


-1 


.8570 


0, 


,1606 


53761, 


,94042 


10, 


,6917 


13, 


,9743 


-1 


.83 


-1 


.5784 


0, 


,1907 


53763, 


,00687 


10, 


,7826 


14, 


,2102 


-1 


.94 


-1 


.6364 


0, 


,1304 


53768, 


,00039 


10, 


,6707 


14, 


,3998 


-1 


.93 


-1 


.7570 


0, 


,1229 


53769, 


,97754 


10, 


,7000 


14, 


,0952 


-1 


.74 


-1 


.6234 


0, 


,3257 


53774, 


,98545 


10, 


,6563 


14, 


,1222 


-1 


.93 


-1 


.6517 


0, 


,2238 


53776, 


,93634 


10, 


,6487 


13, 


,9733 


-2 


.01 


-1 


.5952 


0, 


,1824 


53777, 


,97625 


10, 


.6635 


14, 


,2261 


-2 


.09 


-1 


.6904 


0, 


,1596 


53778. 


,97798 


10, 


.6468 


14. 


,1078 


-1 


.89 


-1 


.6498 


0. 


,1498 


53783, 


,97776 


10, 


.4454 


13, 


,5259 


-1 


.84 


-1 


.4976 


0, 


,1797 


53790, 


,96543 


10, 


,5921 


13, 


,7395 


-1 


.72 


-1 


.5243 


0, 


,1516 


53792, 


,97117 


10, 


,6038 


14, 


,0844 


-1 


.87 


-1 


.6576 


0, 


,2388 


53797, 


,97754 


10, 


.4842 


13, 


,5245 


-1 


.73 


-1 


.4815 


0, 


,2869 


53800, 


,95004 


10, 


.6096 


14. 


,0419 


-1 


.72 


-1 


.6383 


0, 


,1495 



on Hipparcos and Tycho catalogues is presented by the 
VSNET (2005) also in B and V. The UBVRI magnitudes 
were published only for one star "c", which were deter- 
mined by Efimov ct al. (1998) by linking to the standard 
star HD 23949 (Neckel and Chini,'l980). The check star "d" 



Table 5. Characteristics of the prominent minima of the 
"3.2" -hour variability determined using the method of 
"asymptotic parabola" (Marsakova and Andronov, 1996): 
time B JD (and its error) and value of the signal rUmin (and 
its error). 



t, BJD 




at] 


777" 


min 




53728.15648 


0, 


.00180 


10, 


.911 


0, 


.014 


53730.01153 


0, 


.00065 


10, 


.860 


0, 


.009 


53732.00013 





.00122 


10, 


.746 


0, 


.009 


53742.01621 





.00078 


10, 


,854 


0, 


,014 


53742.94917 





.00234 


10, 


,782 


0, 


,019 


53744.00722 





.00220 


10, 


,896 


0, 


,023 


53753.95263 


0, 


.00001 


10, 


,743 


0, 


,007 


53757.00094 


0, 


.00140 


10, 


,846 


0, 


,005 


53757.94716 


0, 


.00072 


10, 


.794 


0, 


.008 


53758.05255 


0, 


.00162 


10, 


.784 


0, 


.010 


53758.97022 


0, 


.00248 


10, 


.747 


0, 


.008 


53760.04089 


0, 


.00030 


10, 


.704 


0, 


.005 



(= GSC 1207 1562), which is closer to TT Ari and suitable 
for measurements, is much redder (B-V=l™15) than "c" 
(B-V=a"69) (Goetz 1985), thus "c" is recommended to be 
used as the comparison star. 

Recently, Henden (2007) published BVRI magnitudes 
for many stars in the field, including "c" (V=ll™05, 
R=10™66) and "d". A compilation of magnitude estimates 
is presented in Table 3. The difference between the bright- 
ness estimates of the comparison star "c" in the filter R 
reaches 0™25. This value is typical for discrepance between 
"standard" stellar magnitudes determined by different au- 
thors. Finally, for the calibration of our observations, we 
adopted the values (V=ll™01, R=10^41) published by 
Efimov et al. (1998). No color correction was made, so 
the instrumental magnitudes are defined as the adopted 
brightness of the comparison star plus instrumental differ- 
ence "variable-comparison". The error estimates for a sin- 
gle measurement of the variable are 0^008 (for the night at 
the top of the outburst) to 0^026 (for 6 nights in the faint 
state). 

3. Results of Analysis of Multi-Component 
Variability 

3.1. Light Curves 

The multi-type character of variability was reported by pre- 
vious researches. In this work, we make a time series anal- 
ysis on an unprecedented number of observational runs, 
which were obtained in the same photometric system R 
(except 1 night in V). The complete light curve, obtained 
during the present observational campaign, is shown in Fig. 
1. According to the character of variability and values of 
the mean brightness m in Table 1, it may be splitted into 
several subintervals: 

- 53668-53700: m ranged from 11^00 to 11^18, and 
this may be characterized as the "pre-outburst" stage. 

However, one may note a slow brightening from 52682 to 
53700 with a best fit slope m = dm./dt = -0.00624(19) 
mag/day. Hereafter the numbers in brackets show the 
error estimate in units of the last digit. Another im- 
portant characteristic of outbursts is a "decay time" 
Tdecay = At / dm for the descending branch of the out- 
burst (cf. Bailey 1975, Warner 1995), and, similarly, a 
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"rise time". Trise = —dt/dm for the ascending branch. 
According to the properties of the least-square linear 
fits, |mTrise| = P^, where p is the correlation coefficient 
between time t and brightness m. For this time interval 
of our observations, Trise = 37.0 ± l.f days/mag. 

- 53700-537f4: "rise to the outburst" , the star brightened 
from m = 11"^09 to 10^26 during U days, T^se = 
16.0 ±0.1 days/mag, m = -0.0584(4). The slope at the 
ascending branch is by 9.4 times larger than in the "prc- 
outburst" stage. The shape of the rise is noticeably lin- 
ear, with a large correlation coefhcient p — —0.9668(65). 
The slope is smaller than that for ourburst of dwarf no- 
vae by an order of magnitude (cf. Warner 1995), thus 
indicating a different nature (possibly an increase of the 
mass transfer rate). 

- 53712-53721: "top of the outburst", two nights 53712, 
53713 from 5 may belong both to the "rise" and the 
"top" , as the brightness continued to increase, whereas 
the value for these two nights fits to the range character- 
istic for the "top" . The mean value is m = 10™ 341 (3). 

- 53721-53723: "upper limit for decay" , During two days 
between these nights, the brightness decreased by 0^"34. 
Because of bad weather in the intermediate night, the 
time interval of 2 days is an upper limit for the total 
duration of decay. 

- 53727-53800: "post outburst slow brightening", the 
mean slope is m = —0.00288(5) mag/day, the values 
of fh for individual nights may differ from the corre- 
sponding linear fit up to » 0^1. 

The observed light curves are shown in Fig. 2. They 

show a variety of shapes and types of variability. Sometimes 
the "3.2— hour variability is prominent (e.g. 53727, 53741- 
53762), sometimes this component of variability has a re- 
markably small amplitude (53685, 53698). The amplitude 
of the quasi-periodic oscillations (QPO), which were also 
previously reported by other authors, varies from night to 
night and also within one run. 

For numerical parametrization of the observed variabil- 
ity, we used several complementary methods of time series 
analysis - periodogram, wavelet analysis, running sine and 
running parabola approximation, 0-C analysis, which are 
described below. 



3.2. Wavelet Analysis 

The wavelet analysis is one of the powerful methods to 

study signals with variable periods and amplitudes (cf. 
Daubechies 1988). For the "'wavelet periodograms" , we 
used the test functions S(P), r{P) (Andronov 1998a) and 
WWZ (Foster 1996), which are specially elaborated for a 
case of finite signals with possibly non-regular arguments. 
Their statistical properties were studied by Andronov 
(1998b). These test functions are complementary, although 
the peaks at them correspond to effective values of the pe- 
riod. Their physical sense may be briefly explained as fol- 
lows (see Andronov (1998ab) for a detailed description): 
r(P) is an effective amplitude (which is equal to a semi- 
amplitude, if the signal is purely sinusoidal); S{P) is a 
weighted mean of the squared correlation coefficient, so 
S{P) « corresponds to a bad periodic fit, and S{P) = 1 
shows that the period is correct (if the analyzed signal is 
purely sinusoidal) ; WWZ is the "energy" signal-to- noise ra- 
tio, so the false alarm "probability" FAPr:! exp(— WWZ). 




Ig P 

Fig. 3. The "wavelet periodograms": test functions used 
S{P), r{P), WWZ for a complete set of observations (see 
Andronov 1998 for a detailed description). The vertical 
lines mark the positions of the highest peaks at S{f), 
which correspond to P = 0^0143 = 20.6 minutes and 
P = 0^135 = 3^*24. 



These test functions are shown in Fig. 3 as a function of 
trial period P. One may note a large scatter of the semi- 
amplitude r(P) at small values IgP (< —2.4 and > 0.0). 
This is an effect of discretncss of times of obscrvions, which 
prevents r{P) to be the best test function among the three 
listed above. The test function S{P) shows highest peaks 
at P = 0^0143 = 20.6 minutes (IgP = -1.844) and 
P = 0^135 = 3^'24 (IgP = -0.87). The corresponding val- 
ues of S{P) (the mean contribution of variability at this 
effective period to the total variance of the signal) are re- 
markably similar: 0.184 and 0.182. The values of WWZ at 
these two best periods are equal to 8.0 and 25.0, respec- 
tively, so both suggested types of variability are statisti- 
cally significant. The mean semi-amplitudes of the period 
are 0^0357 and ffi"0456. The peak at IgP « 0.2 occurs 
at the period, which is comparable with distance between 
the subsequent nights, so the effective number of points is 
small. So we do not interpret this timescale to be of physical 
nature of the object. 

To study stability of the parameters of the variability, 
the "wavelet periodogram" for individual nights are shown 
in Fig. 4. The range of trial periods P is restricted to 3Amin 
to tg — ts, where Amin is the time resolution (time step be- 
tween subsequent observations), and {tg—tg) is the duration 
of observations. Both these parameters are listed in Table 
1. The peaks close to either 20.6— min, or 3.2— hour show 
very significant night-to-night changes both in height and 
the position. 

The period determination of the 3.2— hour variability 
will be discussed in section 3.6, as it exceeds the duration 
of many nightly runs of our observations and consequently 
should not be determined from periodograms for so short 
individual runs because of possible large distortion of the 
smoothing function. 
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Fig. 4. The individual wavelet periodograms S{P) for all nights of observations. To avoid bias, the curves are shown in 
a range of trial periods from 3Amin to tg — tg (see Andronov 1998 for a detailed description). The vertical lines mark 
the positions of the highest peaks of S{f) for the complete set of observations, which correspond to P = 0^0143 = 20.6 
minutes and P = 0^135 = 3^*24. The 5-digit "legend Julian date" LJD is shown near each graph. 
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Fig. 5. Dependence on time of the characteristics of individual runs: sample mean magnitude to, effective magnitude of 
the source of the QPOs niQPo, logarithm of period IgP of QPOs, logarithm of semi-amplitude Igr of QPOs. and the 
corresponding test function S{P). Each point represents one night. 



3.3. Effective Characteristics of QPOs 

We used the following parameters to characterize the QPOs 
for individual nights: the period P, scmi-amplitudc r(P) 
and "effective stellar magnitude" mqpo, which correspond 
to the maximum of S{P) in the interval —2.4 < IgP < 
—1.7. This interval was chosen after the examination of the 
mean and individual "wavelet periodograms" . The param- 
eter TOQPo = m + 2.5 lg((l + 10"-''*'^)/(l - 10-0-4*'-)) is the 
magnitude estimate of the source of QPOs. The physical 
sense is the magnitude of the source of emission, which cor- 
responds to brightening from the stellar magnitude (m -|- r) 
(invisible) to (to — r) (visible). 

These characteristics are listed in Table 4 and shown in 
Fig. 5. The mean semi-amplitudes of variability range from 
(r0139 (run 53760) to 0^0821 (run 53718), i.e. vary by 
a factor of 5.9. The IgP ranges from -2.13 to -1.72, except 
one night 53760 with an outlying value of -2.36 (the run was 
short, so this value represents only a short time interval and 
thus is not very accurate). 

The time variation of to was discussed earlier, with split- 
ting the interval of our observations according to this pa- 
rameter. Surprisingly, the outburst seen in rfi has a very 



different shape and duration in rnqpo- At the beginning of 
the outburst (53703, 53704), when the mean brightness to 
started to increase, the semi-amplitude r and toqpo are 
within their ranges characteristic for the "pre-outburst" 
stage. Then the values of toqpo range from 12™85 to 13™76 
(except one short night 53718) for the nights 53706-53753. 
This stage may be referred to as "QPO-oxcited" , starts 
from the outburst rise and lasts at least 30^^ after the end 
of w 10*^ ourburst. 

After that, rnqpo steeply decreases (i.e. the QPO flux 
increases) with a best fit slope drnqpo / At = —0.0125(37) 
mag/day which is much larger than that of -Afh/dt = 
0.0037(11) mag/day. However, night-to-night variability is 
significant, the slopes only slightly exceed the "3fT— level". 
The best fit for 18 nights (except the short run 53760) is 
TOQPO = 14"05(3) + 2.61(42) • (to - 10"^64), the slope is 
determined at much more significant 6.3ct— level" . This 
argues for a statistically significant relationship between 
mQPo and fh. The slope is larger than unity, so with in- 
creasing mean flux, also increase: the semi-amplitude (re- 
lated to arbitrary amplitude of the flux variations), the 
brightness of the source of QPO (related to the flux of 
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Fig. 6. Dependence on m of the characteristics of individual runs: effective magnitude of the source of the QPOs {rnqpo), 
logarithm of period (IgP) of QPOs, logarithm of semi-amplitude (Igr) of QPOs. Each point represents one night. 



QPO) and the period. The statistical significance of lat- 3.4. QPO Period Variability During the Run 
ter relationship at the 3.4o-— level is seen from the fit 
Igr = -1.628(14) + 0.40(11) • (lgP+ 1.9). 



The relations present above were derived for the part of 
the light curve after the outburst of the QPO. The com- 
plete diagrams showing pairs of parameters (not shown), 
exhibit more complex behaviour and several groups with 
smaller number of points inside. At Fig. 6, the diagrams of 
'rngpo, Ig P and Ig r vs m are shown. One may note two dis- 
tinctly different groups at the intermediate brightness inter- 
val, which are shown by squares (4-th interval) and crosses 
(5-th interval). These time intervals are consequent, and 
are characterized by similar range of ifi. However, they are 
very different by Igr, and, consequently, rnqpo- Contrary 
to the 5-th interval, there are no significant statistical de- 
pendencies on m of any of the parameters rnqpo, ^gP and 
Igr. 



Those relations arc interesting for future modeling of 
accretion disks with slowly increasing brightness. 



Drastic variability of the nightly values of the QPO char- 
acteristics is observed. To study this phenomenon, we have 
used the wavelet test function S(f ) (Andronov, 1998ab), 
as the position of peaks of this fiuiction c;oincides to the 
underlying period. The dependence of the Ig P on time for 
all individual nights is shown in Fig. 7. According to sta- 
bility of the periods during the run, one may distinguish 
between nights with small period variations (e.g. 53712) 
and nndtiple switches (e.g. 53691). Such apparent switches 
at the P{t) diagrams may be explained by aperiodic flares 
or dips, which occasionally change the interval between the 
brightness maxima or minima, respectively. One may ne- 
glect switches at the borders of observations, as the trial 
fit used in the wavelet analysis is filled with observations 
very asymmetrically, what causes larger statistical errors 
and apparent rise or decay of smoothing signal. Examples 
of such jumps of frequency for the "negative superhump" 
state in 1989 were reported by Kraicheva et al. (1999). 
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3.5. Search for Possible Periodic Components 

The presence of QPOs does not mean that there are no 
truly periodic variations, e.g. related to the spin period of 
the white dwarf. Similar situation is present in the mag- 
netic dwarf nova DO Dra, where the new phenomenon - 
"transient periodic oscillations" hides the spin variability 
(Andronov et al., 2008). For this purpose, we have used 
the program "Fo", which realizes the periodogram analysis 
based on the least squares cosine fit (Andronov 1994). The 
test function S{f) is a square of the correlation coefficient 
between the "observed" and "computed" (using the least 
squares sine fit for a trial frequency) values. Statistical sig- 
nificance of the peak increases with its height (see Andronov 
(1994) for description of statistical properties of this test 
function). 

To avoid a bias of the periodogram caused by signifi- 
cant low frequency "3^*2" hour component, it is needed to 

remove it before the periodogram analysis. Similar to BG 
CMi (Kim ct al., 2005) and DO Dra (Andronov ct al., 2008), 
we have used the "running sine" (RS) fit, with an initial 
period P = 0^01413 taken from the wavelet periodogram. 
Generally, the analysis of the characteristics of this fit (not 
shown) confirms our previous findings - variability of the 
period and semi-amplitude. Thus we have used only esti- 
mate of the low-frequency trend to compute the residuals 
of the signal from the smoothing values (0-C). 

The test function S{f) for the residuals (0-C) is plotted 
versus frequency in cycle/day. For reference, the frequency, 
which corresponds to the best period corresponding to the 
mean wavelet periodogram, is 69.8 cycles/day. The indi- 
vidual periodograms are shown in Fig. 8. At some nights, 
the peaks are high and narrow (e.g. 53692, 53731, 53769), 
which correspond to oscillations, which are stable during 
the night. Wide peaks are often high (53700, 53707, 53718) 
and wide, what is consistent with the inverse proportional- 
ity of the peak width to the duration of observations. The 
periodograms show multiple peaks (e.g. 53668), what cor- 
responds to absence of stability of the periods during the 
nights. The frequency corresponding to the highest peak, 
is a good parameter for QPOs. However, due to period 
changes, the highest peak may not correspond to that at 
the "wavelet periodogram" S{P), and the estimate of the 
semi-amplitude r is much smaller than that for the wavelet 
(they coincide only, if the signal is strictly periodic). Thus 
we continue to use estimates from the wavelet analysis. The 
periodogram for all nights (not shown because of a huge 
length of the file) shows the highest peak at the frequency 
/ = 68.2712(3) (period P = 0'?0146475(6) = 21.0924 min- 
utes, semi-amplitude r = 0^^107(6), initial epoch for the 
maximum brightness Tmax = 53731.73880(14). Although 
S{f) = 0.02 is small, i.e. only two per cent of variability 
may be explained by this periodic component, the "false 
alarm probability" FAP=10~^^'^ is very small, indicating 
that this variability is statistically significant. 

The estimate of the semi-amplitude is much smaller 
than that from the wavelet analysis, so this peak may be 
resulted from several nights with QPOs with similar fre- 
quency. 

3.6. Superhump Period 

Determination of the period is complicated, as the sys- 
tem exhibits several types of variability. The periodogram 



analysis needs a preliminary removal of the slow trend, 
which is not well defined because of the outburst and 
mean brightness fluctuations. Moreover, the light curve 
sometimes shows two maxima per preliminary siiggested 
"3l'2-hour" period. Thus we have used the times of promi- 
nent minima of the light curve, following an approach 
used during our previous campaigns (Tremko et al., 1996, 
Andronov et al., 1999). For better statistical accuracy, we 
used the method of "asymptotic parabola" (Marsakova and 
Andronov, 1996). The list of 12 minima is presented in 
Table 5. 

To seach for possible periods, we analyzed these 12 times 
of minima. The program "PerMin" was used, which realizes 
the "method of characteristic events" (Andronov, 1991). In 
this method, two parameters (period and initial epoch) are 
being determined, which minimize the r.m.s. deviation of 
phases from zero. The corresponding dependence of the test 
function cr{P) on the trial period is shown in Fig. 9. 

The best fit elements are: 

Tmin = BJD2453747.0700(47) -|- 0.13232(5) • E (1) 

The period value is close to that published for previous 
campaigns: P = 0':'l32959(13) (Tremko et al., 1996) and 
0"? 133160(4). The difference between the pairs exceeds the 
statistical errors, but all these values are much smaller, 
than that in the "positive superhump" state P = 0^14926 
(Skillman et al., 1998) or 0^14840 (Andronov et al. 2005). 
So the star during our set of observations was in the "neg- 
ative superhump" state. 

The deviations of phases from zero sometimes reach 0.2, 
what prevents the interpretation of these minima as peri- 
odic eclipses. We tried to find a better ephemeris by remov- 
ing some of the minima, but it has not resulted in better 
coincidence of phase curves for all nights. 

The phase light curves smoothed using the method of 
"running parabola" (Andronov, 1997) are shown in Fig. 10. 
The time is expressed as a continuously increasing phase 
computed using the ephemeris (1), i.e. from 1 to 3, not as 
usually defined from to 1. So practically the brightness 
variability is shown vs. time (divided by the period) and 
shifted by some integer number. In TT Ari, the instability 
of the phase light curve is present not only from night to 
night, but also from one cycle of variability to another one. 
So this kind of presentation of the light curve is preferred, 
like that used for the intermediate polars (e.g. Kim et al., 
2005). 

The filter half width is At = 0'?05, the value adopted 
by Tremko et al. (1996) and Andronov et al. (1999). The 
positions of many minima, which were not so prominent, 
and thus not listed in the Table 5, are in a good agreement 
with the ephemeris above. However, for other nights, at 
phase zero, there may be observed a secondary minimum 
or even a maximum. The light curves may be classified as 
"one-maximum" or "two-maxima" per period. 

This result may be explained by strong variations of the 
disk luminosity, possibly due to the inhomogeneity. 

4. Conclusions 

Based on the analysis of 51 observational runs (JD 
2453669-2453800), which were obtained at the 40-cm tele- 
scope of the Chungbuk National University (Korea) and 
cover 226 hours, the following conclusions are made: 
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Fig. 7. (Available electronically only). Dependence of the local value of the period P of the quasiperiodic oscillations 
(QPO), which corresponds to a highest maximum of the test function WWZ in the adopted interval —2.4 < IgP < —1.7, 
for the individual nights of observations. The short horizontal lines mark the value of IgP = —1.844, which corresponds 
to the maximum of the wavelet periodogram S{P) for all data. The ordinate is marked from -2.4 to -1.7 for all nights. 
The 5-digit "legend Julian date" LJD is shown near each graph. Time is expressed as BJD-LJD, like in Fig. 2. 
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Fig. 8. The pcriodograms S{f) for the individual runs. The 3.2-hour cychcal trend was removed using the "running sine" 
method. The 5-digit "legend Julian date" LJD is shown near each graph. 
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Fig. 10. (Available electronically only). Fits for individual nights using the method of running parabola with the filter 
half- width At = (fi05. The time is expressed as a continuously increasing phase computed using the ephemeris (1). The 
5-digit "legend Julian date" LJD is shown near each graph. 
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Fig. 9. The dependence of the test function <j{P) for 12 
selected moments of minima on the trial period P. The 
position of the lowest minimum at P=0'? 132322 is shown 
by the short line. 



tometric behaviour of the system is dependent on the 
current mean magnitude, but also on the previous state. 

— The analysis of the pairs of characteristics separately 
for these groups shows statistically significant correla- 
tions only for the state "5" . For this state, with increas- 
ing mean brightness, also increase the period, amplitude 
and brightness of the source of QPOs. 

— The brightness variations during the present stage of 
"negative superhumps" range from 10™27 to 11™18, 
whereas the mean brightness in the "positive super- 
hump" state was only slightly smaller: R~ 11™3 
(Andronov et al. 2004). This means that even small 
variations in mean brightness (few tenth of magnitude) 
are sufficient to switch between "positive" and "nega- 
tive" superhumps. "Positive" superhumps are excited 
at lower luminosity, than the "negative" ones. As TT 
Ari sometimes exhibits low states to « 17™, it is an 
interesting observational task to study, which type of 
variability dominates during low luminosity state. 
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— During the period of observations, TT Ari was in the 
"negative superhump" state. The preceeding "positive 
superhump" lasted from 1996 to 2004. 

— As in the previous "negative superhump" states, the 
light curve is characterized by the "3.2-" hour vari- 
ability and a superimposed quasi-periodic oscillations. 
The shape of the "superhump" curve is highly variable, 
showing one-maximum (per period) or two-maxima 
structure. 

— The ephemeris for 12 best pronounced "super- 
hump" minima is T„„ = ^7132453747.0700(47) + 
0.132322(53)i?. This is in the range detected in previous 
studies for the "negative superhumps" . The wide range 
of the phases of minima (from -0.2 to 0.2) indicates that 
these minima may not be eclipses of the source of emis- 
sion, which is fixed in the rotating frame, but the po- 
sition may be dynamically variable. Another source of 
the distortion of the "superhump" light curve (and, con- 
sequently, of the phase scatter) is a strong flickering or 
aperiodic small- amplitude flares. 

— The quasi-periodic oscillations (QPO) are present with 
a mean "period" of 21.6 min and mean semi-amplitude 
of 36 mmag. The position and height of the peaks at 
the "wavelet periodogram" vary from night to night and 
sometimes during the night. 

— The dependence on either time, or the mean bright- 
ness, of the magnitude of the source of QPOs, pe- 
riod and amplitude shows that the interval of observa- 
tions was splitted into 5 parts, showing different char- 
acteristics, namely: 1) the "pre-outburst" stage, 2) the 
"rise to outburst" , 3) "top of the outbursts" , 4) "post- 
outburst QPO" state, 5) "slow brightening". The new 
phenomenon was detected, which we called the "post- 
outburst QPO" state. This state is characterized by a 
large amplitude of the QPOs for more 30 days ("4") 
after the outburst (part "3") lasting 10 days. During 
both states "3" and "4" , the amplitude of the QPOs 
was larger, than in the states "1" and "5" . So the pho- 
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